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We examine the changes in kinetic energy dissipation of a turbulent channel flow caused 
by a spanwise magnetic field. The numerical study is based on our simulation data from 
[I] obtained by direct and large eddy simulations. We find that the Joule dissipation 
can exceed the viscous dissipation in the weakly dissipative bulk region, but remains 
comparatively small in the turbulence-generating near-wall region. 

Introduction This paper continues our recent analysis [11 of turbulent 
channel flow in a steady magnetic field imposed in the spanwise (parallel to the 
walls but normal to the mean flow) direction. The magnetic Reynolds number 
is assumed small, and the problem is solved in the framework of the quasi-static 
approximation j2j. We are interested in this problem because the channel flow 
with spanwise magnetic field represents the simplest natural flow configuration, in 
which we can study turbulence under the combination of two conditions common 
in industrial and laboratory flows of liquid metals: imposed magnetic field and 
mean shear. At the same time, the transformation of turbulent fluctuations oc- 
curs without direct interaction between the mean flow and the magnetic field and 
without the dominating influence of Hartmann boundary layers. 

A thorough study of the flow at high hydrodynamic Reynolds number Re = 
UL/v = and moderate Hartmann number Ha = BL{a / pv)^^^ was conducted in 
[1]. Here, U and L are the typical velocity and length scales, chosen in our case 
as the centerline velocity of the laminar flow and half the channel width, B is the 
strength of the applied magnetic field, and cr, p and v are the electric conductivity, 
density, and kinematic viscosity of the fluid. The results confirmed and extended 
the conclusions of the earlier experimental work and numerical simulations 
at low Reynolds number [4 . We found that the main effect of the magnetic 
field is the suppression of turbulent fluctuations, particularly pronounced in the 
wall buffer regions, where the fluctuations are generated by strong mean shear. 
The resulting reduction of the turbulent momentum transport in the wall-normal 
direction decreases the wall friction drag and transforms the mean flow profile, 
which becomes steeper, acquires higher centerline velocity, and loses the zone of 
the logarithmic layer behavior. 

Analyzing the anisotropy of turbulent fluctuations, we found that the elonga- 
tion of flow structures in the direction of the magnetic field observed in homoge- 
neous zero-shear turbulence (see, e.g. [5j|6l[7]) appears only in a weak form and 
only in the central area of the channel. It should be stressed that all the results 
of [T] were obtained at moderate strength of the magnetic field, at Ha < 30 for 
Re = 10^ and Ha < 40 for Re = 2x 10"*. Stronger magnetic fields lead to an inter- 
esting intermittent flow regime, in which turbulent states alternate with periods 
of nearly laminar behavior . 



The combination of the two factors, the mean shear and the imposed magnetic 
field, also renders the channel flow an excellent configuration for testing computa- 
tional models of turbulence, in particular LES models. This was done in [T] using 
the a-posteriori comparison with the results of high-resolution DNS. In agreement 
with the earlier studies, such as [Hlinill^i the dynamic Smagorinsky model was 
shown to accurately reproduce the MHD flow transformation. 

In the present paper, we extend the investigation by considering the viscous 
and magnetic (Joule) dissipations. The databases generated in the LES and high- 
resolution DNS computations 1] are employed to analyze the effect of the magnetic 
field on the dissipation rates in various zones of the flow. An important goal of the 
study is to interpret the observed flow transformation in terms of the local energy 
balance. We provide further verification of the dynamic Smagorinsky model by 
analyzing the accuracy, with which the model reproduces the dissipation rates. 
Such accuracy is generally considered one of the most important characteristics of 
an LES model's performance. 



Equations and numerical method We consider the flow of an incom- 
pressible, electrically conducting fluid in a plane channel between two insulating 
walls located at z = ±L. The flow is driven by a streamwise pressure gradient and 
submitted to a constant spanwise magnetic fleld B — Bey . The mean flow velocity 
Uq is kept constant in the simulations. The dimensional governing equations and 
boundary conditions for the velocity field u.^, pressure p, electric current density 
Ji and electric potential are 
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In these equations, Ty — 2pi'Sij is the resolved viscous stress tensor, where 
2Sij = dui/dxj + duj/dxi is the rate-of-strain tensor. The additional subgrid 
stress tensor Tj", which appears in LES computations, is given by r° = 2pvTSij. 
It is computed using the local filtered rate of strain tensor Sij, with the eddy 
viscosity vt determined by the dynamic Smagorinsky model with averaging in 
horizontal planes. 

Variables are made non-dimensional by the laminar centerline velocity [/, half- 
channel width L and ULB (for 0). The non-dimensional parameters are Re and 
Ha. The problem is solved numerically in a domain with the periodicity lengths 
Lx/L — 27r and Ly/L = tt. The pseudo-spectral numerical method is described 
in [U [TT] . The numerical resolution is x Ny x — 256'^ collocation points in 
the simulations with Re = 10'' and 512^ x 256 collocation points at Re = 2 x 10'*. 
The spatially filtered equations for LES are solved using 64'^ and 128'^ collocation 
points at Re ~ 10'' and Re = 2 x 10'', respectively. 

The equation of the local kinetic energy balance is formed from the momentum 
equation ([1]) by multiplying it by the velocity vector tt, making use of Ohm's law, 
and rearranging such that the energy transport terms are entirely separated from 
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the local dissipation terms. The result is the equation 



de 
di 



_ = -V • Af - e^igc - - esgs, (5) 



where e = is the specific kinetic energy, 

^visc = '^ij ^ij I P (^) 
is the viscous dissipation rate, and 

= J^I<^P (7) 

is the Joule dissipation rate. The additional dissipation rate due to the subgrid 
stresses is 

esgs = TtfS,,lp. (8) 

A/' is the energy flux vector due to convective transport, pressure force, friction 
force, subgrid stresses, and Lorentz force. It is given by 

pN^ = pu^ ( Y ^ p ) ^ "'^'^ ^ "'^5 ~ 

In the case of DNS, the subgrid stress contribution is absent from all of the above 
equations. 



Results We first consider the spatial structure of the Joule dissipation field. 
Fig. [1] illustrates the trend observed in our earlier paper [T] for the velocity field: 
a stronger magnetic field leads to turbulence suppression and smoother, larger 
structures. 
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Figure 1: Iso-surfaces E/^/e™* = 2 of the Joule dissipation rate for Ha = 10 
(left) and Ha = 30 (right) from DNS with Re 10000. Snapshots of the entire 
computational domain are shown. 



For further analysis, we focus on the average distributions of dissipation rates 
and their transformation by the magnetic field. The averaging is performed in 
a;-y-planes and in time over several convective times L/U. We begin with the 
demonstration that LES reproduces the dissipation rates with satisfactory accu- 
racy. Fig. [2] shows close agreement between the mean profiles of the viscous 
dissipation rate from DNS and the sum of viscous and sub-grid scale dissipation 
rates from LES in the non-magnetic cases. Figs. [3Ia,b) show similar agreement 
for the cases of the strongest magnetic fields realized in the simulations. For the 
latter two cases, the Joule dissipation rates shown in Figs. [3jc,d) also show good 
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Figure 2: Time-averaged profiles of Eyisc ^'^^ (a) Re = 10000 at Ha = and (h) 
Re — 20000 at Ha = 0. In case of LES tfie total viscous dissipation ^visc+sgs 
plotted. Normalization is by u^/v obtained for each case from DNS. The distance 
z"*" from the wall is measured in units of the friction length v/ut. 
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Figure 3: Time-averaged profiles of Eyisc '^'^^ ^t^ (^''^) ~ 10000 at Ha = 30 
and (h,d) Re = 20000 at Ha = 40. In case of LES the total viscous dissipation 
Eyisc-l-sgs ^'^'^ resolved part of the Joule dissipation gresolved plotted. 
Normalization is by u^/v obtained for each case from DNS. 



agreement. It should be mentioned that the comparison should, in principle, be 
done between the LES and filtered DNS results. Filtering was not performed in 
our analysis. Unlike the viscous dissipation, the Joule dissipation rate varies with 
the length scale in about the same way as the kinetic energy, and the Joule dis- 
sipation in the subgrid scales should therefore be quite small. The slightly larger 
discrepancies in near the wall at Re = 10^ in Fig. [3Jc) are most likely caused 
by a comparatively coarser LES resolution than for Re — 2 x (Fig- Eld)). 
In any case, our LES results are quite satisfactory for the mean profiles of the 
dissipation rates. This observation is further supported by Table [U which shows 
that the fraction of the Joule dissipation relative to the total dissipation is well 
represented in the LES model. 

Non-magnetic turbulent channel flows are characterized by local equilibrium 
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Table 1: Time- and volume-averaged subgrid-scale and Joule dissipation rates 
relative to the total energy dissipation rate e^Q^ = eyisc + ^/j + ^sgs- 



between kinetic energy production and dissipation and by the approximate relation 
^visc ~ 1/2^ in the logarithmic layer !12j. The product ^^^Eyisc should therefore 
be approximately constant. Fig. [3] shows this quantity as found in our DNS 
computations. We can identify a plateau for Re — 2 x in both magnetic 
and non-magnetic cases, although it is shown in 1 that the mean velocity profile 
for Ha = 40 deviates significantly from the logarithmic law. For Re = 10'* the 
quantity Eyisc decays somewhat faster than l/z+ in the non- magnetic as well as 
the magnetic case, and the plateau is therefore not observed. 




Figure 4: Compensated viscous dissipation z^SyigQ from DNS for Re ~ 10000 at 
Ha = 0,30 and Re = 20000 at Ha = 0,40. The plots are continued up to the 
middle of the channel, which varies because of different friction lengths v /ur- 



The systematic trends in the behavior of the dissipation rates with increasing 
magnetic field are illustrated in Fig. [5] by the LES data obtained at Re = 10**. 
The results for Re = 2 x 10'* are qualitatively identical and, therefore, omitted. 
Figs. [3a,b) show that the viscous dissipation rate is significantly reduced by the 
magnetic field near the wall and in the middle of the channel. The reduction at 
the walls is expected because of the suppression of the intensity of the turbulent 
fluctuations observed in [1] . The reduction in the middle is less signiflcant for the 
overall level of dissipation since this region contributes very little, but the relative 
change in £yjgQ by a factor of about 3.3 from Ha = to Ha — 30 is larger than at 
the wall. The Joule dissipation rate shown in Fig. [DJc) increases significantly with 
Ha but without much change in the profile shape. The subgrid-scale dissipation 
rate in Fig. [DJd) decreases significantly from Ha = to Ha = 30, as can be 
expected from the magnetic damping of the turbulence and the corresponding 
depletion of the subgrid scales. 
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Figure 5: Horizontally and time-averaged rates of (a,b) viscous, (c) Joule and 
(d) subgrid-scale dissipation rates from LES for different Ha at Re = 10000. The 
viscous term is shown for the near- wall region (a) and, to facilitate the comparison 
in the bulk region, in log-scale (b). Quantities are normalized by U^/L, which is 
identical for the different Ha. 



The relative contributions of the dissipation rates are compared in Fig. [SI 
which shows mean profiles calculated in LES for Re — 10* and Re = 2 x 10*. 
At low Ha, the contribution from the Joule dissipation rate is everywhere smaller 
than that of the viscous dissipation, but for the highest Ha, the Joule dissipation 
becomes comparable or even larger than Eyisc bulk region. The overall 

contributions are compared in Table [TJ The contribution of the subgrid-scale 
dissipation diminishes considerably in the bulk with increasing Ha, especially for 
Re — 2 X 10* with its larger number of grid points. The overall contribution of 
the subgrid-scale dissipation remains nonetheless of the order of 10%. 

In conclusion we remark that both the Joule and viscous dissipation rates 
peak at the wall. The Joule dissipation rate attains a local maximum in the 
buffer region, where the turbulence is primarily generated. The shape of the Joule 
dissipation profile is not much affected by the strength of the magnetic field. Its 
contribution to the total dissipation remains fairly modest, but it can exceed the 
viscous dissipation locally (in the bulk) before the intermittent dynamics described 
in [5] appears at larger Ha. The contribution of the subgrid-scale dissipation 
and, thus, the importance of the LES modeling, diminishes with the strength of 
the magnetic field, but remains significant even at the strongest field at which a 
continuous turbulent state is maintained. 
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Figure 6: Mean profiles of the dissipation rates obtained from LES for Re = 10000 
(left) and Re = 20000 (right) at different Ha. Normalization is hy ul^/v obtained 
for each case from DNS. 
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